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Abstract 

Half-metallic ferromagnets like the full Heusler compounds with formula X2YZ are supposed 
to show an integer value of the spin magnetic moment. Calculations reveal in certain cases of 
X = Co based compounds non-integer values, in contrast to experiments. In order to explain 
deviations of the magnetic moment calculated for such compounds, the dependency of the electronic 
structure on the lattice parameter was studied theoretically. In local density approximation (LDA) , 
the minimum total energy of Co2FeSi is found for the experimental lattice parameter, but the 
calculated magnetic moment is about 12% too low. Half-metallic ferromagnetism and a magnetic 
moment equal to the experimental value of 6ub are found, however, only after increasing the lattice 
parameter by more than 6%. 

To overcome this discrepancy, the LDA+U scheme was used to respect on-site electron correlation 
in the calculations. Those calculations revealed for Co2FeSi that an effective Coulomb-exchange 
interaction U e ff = U — J in the range of about 2eV to 5eV leads to half-metallic ferromagnetism 
and the measured, integer magnetic moment at the measured lattice parameter. Finally, it is shown 
in the case of Co2MnSi that correlation may also serve to destroy the half-metallic behavior if it 
becomes too strong (for Co2MnSi above 2eV and for Co2FeSi above 5eV). These findings indicate 
that on-site correlation may play an important role in the description of Heusler compounds with 
localized moments. 

PACS numbers: 71.20.Lp, 75.30.Cr, 75.50.Cc 
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I. INTRODUCTION 



Half-metallic ferromagnetism was first predicted for the half Heusler compound NiMnSb 
by de Groot et a^. A half-metallic ferromagnet is a material with a gap between spin states 
at the Fermi energy (e^), such that the majority states are conducting and the minority 
states are semi-conducting. Half-metallic ferromagnets (HMF) are supposed to exhibit a 
real gap in the minority (or majority) density of states (DOS). This gap has in compounds 
with integer site occupancies the consequence that the magnetic moment has to be an integer, 
because the number of occupied minority states is integer. It should be noted that the HMF 
character is lost if an already small deviation from an integer value is observed. 

Certain X 2 YZ Heusler compounds have also been predicted to show half-metallic 
ferromagnetism 2 . Co and Mn containing Heusler compounds have attracted particular at- 
tention, as they are strongly ferromagnetic with high Curie temperatures (above 900K)*- 
One of them - Co2MnSi - is an attractive candidate for many applications in spin dependent 
electronics. This compound was used by several groups to produce thin film o 4l5iei l8 > 9 and 
device o 10 ^ 1 . It was recently shown that Co2FeSi is also a very good candidate as it exhibits 
a very high Curie temperature of HOOK and a magnetic moment of 6/15 12 . This compound 
is somewhat peculiar regarding its Curie temperature and magnetic moment because both 
are reaching the highest values observed up to now in this class of materials. 

In Heusler as well as other compounds, the half-metallic ferromagnetism has not yet been 
unambiguously confirmed by experiments, although there is strong evidence in some cases. 
Despite a lot of work on Co and Mn containing Heusler compound a 5 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 ' 19 ' 20 ' 21 ' 22 , a 
breakthrough is still required to prove the existence of half metallic ferromagnetism in these 
materials. 

The electronic structure plays an important role in determining the magnetic properties of 
the Heusler compounds. In this work, a comprehensive investigation of the equilibrium struc- 
tural, electronic and magnetic properties of Co2FeSi is presented and compared to Co2MnSi. 
The dependence of the magnetic moment on the lattice parameter was analyzed, focusing 
on the differences between LSDA (local spin density approximation) and GGA (generalized 
gradient approximation) treatments. A series of calculations for C02 based Heusler com- 
pounds was performed which shows that mostly all C02 based Heusler compounds exhibit 
half-metallic ferromagnetism 23 , in agreement with Ref.— . The Slater-Pauling rul o 12124125126127 
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fits for most C02 based Heusler compounds and may be used as a rule of thumb for esti- 
mating their magnetic moments. According to this rule the magnetic moment of Co2FeSi is 
expected to be 6/x B per unit cell. 

This work was stimulated by a recent experimental study that was devoted to the mag- 
netic moment and Curie temperature of the Co2FeSi compound 12 . To provide a concise 
picture of the Co2YSi systems (Y = Mn, Fe), a theoretical investigation will be presented 
utilizing the local (spin) density approximation (L(S)DA), GGA, and LDA+?7 methods. 

II. COMPUTATIONAL DETAILS 

Heusler compounds 28 ' belong to a group of ternary intermetallics with the stoichiometric 
composition X2YZ ordered in an L2i-type structure, space group Fm3m, many of which are 
ferromagnetic 29 . Remarkably, the prototype CU2M11AI is a ferromagnet even though none 
of its constituents is one^. In general, the X and Y atoms are transition metals and Z is a 
main group element. In some cases Y is replaced by a rare earth element. The X atoms are 
placed on 8a (1/4,1/4,1/4) Wyckoff positions and the Y and Z atoms on 4a (0,0,0) and 4b 
(1/2,1/2,1/2). The cubic L1\ structure consists overall of four interpenetrating fee lattices, 
two of them equally occupied by X. The latter two X 2 fee sub-lattices combine to a simple 
cubic sub-lattice. Thus, the, compound may be seen as a CsCl (B2) superstructure with 
the cubic vacancies of the X2 sub-lattice filled alternatingly by Y or Z atoms. Other than in 
a simple CsCl structure, the nearest neighbor polyhedra of each X atom are given by two 
different tetrahedra (rotated by 90° with respect to each other) that are built from either 
Y or Z atoms. The T point of the paramagnetic structure has the symmetry Oh- However, 
the wave functions at the T point have to be described by m the ferromagnetic state to 
account for the correct transformation of the electron spin. 

Self consistent band structure calculations were carried out using the full potential linear 
augmented plane wave (FLAPW) method as provided by P. Blaha et al (Wien2k) 31 . The 
calculations were performed using the von Barth-Hedin^ 2 - parametrization of the exchange - 
correlation functional within LDA. Further, the Perdew-Burke-Ernzerhof^ implementation 
of GGA was used. GGA accounts specifically for density gradients that are neglected in 
pure L(S)DA. The energy threshold between core and valence states was set to -6Ry. We 
considered 2.?>a.B hr for the muffin-tin radii (Rmt) of all atoms leading to nearly touching 
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spheres. 455 irreducible k points of a 25 x 25 x 25 mesh, were used for Brillouin zone 
integration. The number of plane waves was restricted by Rmt x k max = 7. The energy- 
convergence criterion was set to 10 -5 . Simultaneously, the charge convergence was monitored 
and the calculation was restarted if it was still above 10~ 3 . 

Finally, the LDA+U method^ was used to account for on-site correlation at the transition 
metal sites. The LDA+U method accounts for an orbital dependence of the Coulomb and 
exchange interaction that is absent in pure LDA. Here, the effective Coulomb-exchange 
interaction (£7 e // = U — J) was used for the calculations. This particular scheme is used 
in Wien2k to include double- counting correction. One leaves the framework of density 
functional theory (DFT) by applying the LDA+U method. However, it is one of the most 
popular and easiest approaches to consider electron-electron correlation. It will be shown 
that the LDA+U method gives qualitative and quantitative improvements compared to the 
bare LDA or GGA approaches, when applied to Co2FeSi. 

III. RESULTS AND DISCUSSION 

The main focus was put on the magnetic moment as this shows a very strong discrepancy 
if comparing the experimental and calculated values reported for Co2FeSi. Niculescu et ati& 
reported already a magne, tic moment of 5.9/zb per unit cell at 10K, whereas band structure 
calculations using the screened Korringa-Kohn-Rostocker method (KKR) predicted later 
only 5.27/iH^i. The present LSDA calculations (FLAPW) revealed also a total magnetic 
moment of only 5.29 \xb- This value is much lower than the measured value of about 
Such a large discrepancy points clearly that the actual electronic structure of this compound 
is different from the calculated one. Even so the calculations for Co2MnSi agree with the 
experiment concerning the magnetic moment, their reliability concerning the details of the 
electronic structure may have to be called in question as long as the differences obtained for 
Co2FeSi stay unexplained. 

The discrepancy between experiment and calculation was not removed when using the 
GGA parametrization resulting still in a too low value (5.56/ie). Inspecting the spin resolved 
density of states (DOS) and band structure revealed the appearance of a gap in the minority 
bands, but located below the Fermi energy (see FigslHbEfe i n Sec Jill B"2|) . This is finally one 
reason why the magnetic moment is too low and not an integer. In the following, details of 
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the electronic structure calculations will be discussed. 
A. Structure optimization 

First, a structural optimization was performed for Co2FeSi and Co2MnSi to determine 
whether the experimental lattice parameter minimizes the total energy. It was found that 
the optimized lattice parameter from the calculation agrees very well with experimental 
values of a exP) Fe = 5.64Aand a exP) Mn = 5.645A. Within the FLAPW scheme the structure 
optimization for Co2YSi (Y = Mn, Fe) was performed using the GGA parametrization and a 
spin polarized (ferromagnetic) setup. The initial crystal structural parameter for Aa/a = 
is taken from the experimental value a exp and varied in steps of 2%. For each of the different 
a parameter, the iteration procedure was performed up to self-consistency. The energy 
minimum then defines the optimal value of a. 

For Co2FeSi, the energy minimum (ferromagnetic) was found to appear at a = 
5.63A(corresponding to Aa/a exp « —0.2%, see FigO] GGA). It was also found that the 
ferromagnetic (spin polarized) configuration has the lower minimum total energy compared 
to the paramagnetic (non spin polarized) case, not shown here. Similar calculations re- 
vealed for Co2MnSi that the energy minimum appears at a = 5.65lA(corresponding to 
Aa/a exp « +0.1%. 

Figure [T] displays in addition the results of a structural optimization using the LDA+U 
scheme with U e ff as variational parameter. The results of this optimization being performed 
for different values of U e ff will be discussed below in Sec JIII Bl 

1. Electronic structure and lattice parameter 

In order to understand why LSDA and GGA do not give the expected magnetic moment 
and position of the gap, the dependence of both on the lattice parameter was carefully 
inspected for Co2FeSi as well as Co2MnSi. It is important to compare the electronic structure 
of Co2FeSi with Co2MnSi as the latter is predicted to be a half-metallic ferromagnet with 
the measured magnetic moment of 5/i# found in calculations using the experimental lattice 
parameter. However, a spin polarization of only about 55% was obtained for this compound 
as determined by point contact Andreev reflection spectroscopy^^- The spin polarization 



5 




FIG. 1: Structural optimization for Co2FeSi. 

Shown is the change of the total energy as function of the lattice parameter (ao = 5.64A). The 
GGA results are compared to LDA+U calculations with different U e jf (see text). Note that all 
energy scales are shifted to result in Et o t(0) = for better comparison. The lines are results from 
a polynomial fit. 

at the Fermi energy determined by means of photo emission was even lower (8 . . . 11 %) 8 . 

The calculated total and site specific magnetic moments are shown in Fig|2]for Co2FeSi 
and Co2MnSi. In Co2FeSi (Figl^)), the calculated atomic resolved magnetic moments of 
Co and Fe increase both with a and the overall magnetic moment follows the same trend, 
it increases from 4.96/is to above 6/ib with increase of a. The experimental value of 6fiB 
is found for an enlargement of the lattice parameter by about 6 . . . 10%. The site specific 
moments seem to saturate above +6% change of the lattice parameter at about 1.5/ig and 
3/xb for Co and Fe, respectively. (Slightly higher values are compensated by an anti-parallel 
alignment of the moments at Si sites and in the interstitial.) 

The total magnetic moment of Co2MnSi (Fig(2K) increases also slightly with the lattice 
parameter but it stays at 5/xb in the range of ±6% change of a. There is a significant 
change in the magnetic moment at Mn sites, it increases with lattice parameter. At the 
same time, the magnetic moment of Co decreases. Thus the Co moment counter balances 
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a) Co 2 MnSi b) Co^FeSi 
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FIG. 2: Lattice parameter dependence of the magnetic moments. 

Shown are the total and site specific magnetic moments of Co2MnSi (a) and Co2FeSi (b) as function 
of the lattice parameter. Experimental values are assigned by a cross. The experimental lattice 
parameters are ay=Mn = 5.645Aand ay=Fe = 5.64A. Lines are drawn through the calculated 
values to guide the eye. 

the Mn moment such that the overall magnetic moment keeps constant while changing the 
lattice parameter. 

At this point one can say that a moderate change of the lattice parameter does not change 
the overall magnetic moment of Co2MnSi. The case of Co2FeSi is completely different and an 
overall change of about 1/zb is observed in the same range of Aa/a as for Co2MnSi. Overall, 
the magnetic moment of Cc^FeSi is less stable against variation of the lattice parameter 
(at least if using the same parameters for integration and convergence criteria as for the 
Co2MnSi calculations). It exhibits some fluctuations about the integer value at very large 
lattice parameter. 

The formation of the gap and localized magnetic moments in Heusler compounds is due 
to hybridization as was explained in detail by Kiibler et a f 7 ^ 38 . There is a close relation 
between the magnetic moment and the HMF character. The appearance of the gap in the 
minority density constrains the number of minority electrons to be integer—. However, an 
integer value of the magnetic moment may not result automatically in a real gap in the 
minority (or majority) density. The band structure of Cc^FeSi (see also below in Sec lIII B"2|) 
revealed already a small gap in the minority states in the calculation for the experimental 
lattice parameter, but this gap is located below the Fermi energy (ep)- Therefore, the band 
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FIG. 3: Dependence of the minority band gap on the lattice parameter. 

Shown are the extremal energies of the gap involving states for Co2MnSi (a) and Co2FeSi (b). The 
shaded areas assign the region of half-metallic ferromagnetism. Lines are drawn to guide the eye. 

structure was closer inspected in order to proof whether the integer moment is related with 
the appearance of a real HMF minority gap. 

Figure H3 shows the dependence of the extremal energies of the lower (valence) band and 
the upper (conduction) band of the minority states enveloping the gap. For both materials 
the magnetic moment has to be integer in the region were ep falls into the gap (gray shaded 
areas in FigOJ), that is the region of half-metallic ferromagnetism. 

Figure shows the dependence of the gap on the lattice parameter for Co2MnSi. The 
shaded area corresponds to the region of half-metallic ferromagnetism. It is clearly seen that 
Co2MnSi behaves like a HMF within ±5% change of the lattice parameter. Therefore, a 
moderate change of the lattice parameter does not change the HMF behavior of Co2MnSi. 
Figure shows the dependence of the gap on the lattice parameter for Co2FeSi. It is 
seen that the gap encloses the Fermi energy between 7.5% and 12% enlargement of the 
lattice parameter. The gap is completely closed for lattice parameter being 4% smaller 
than a exp . Obviously the gap (AE max = 0.4eV at +12%) is smaller compared to C^MnSi 
(AE max = 0.9eV at -2.5%) and appears in a narrower range of the lattice parameter. So 
far it is seen that the HMF character and thus the integer magnetic moment is more stable 
against variation of the lattice parameter in Co2MnSi compared to C^FeSi. 

An enlargement of the lattice parameter by 8 . . . 12%, as needed to explain the magnetic 
moment of Co 2 FeSi, corresponds to a volume-expansion of about 26 . . . 40%. Such a large 
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FIG. 4: Electronic structure of Co2FeSi with enlarged lattice parameter. (LSDA-GGA, a = 6.2A) 

expansion of the crystal volume by about 1/3 is rather unrealistic and falls far out of the ex- 
pected uncertainties for the experimental determination of a. However, it may be interesting 
to inspect which changes in the electronic structure are caused by such an expansion. 

The density of states (DOS) and the band structure of Co2FeSi are shown in Fig0]for the 
enlarged lattice parameter. The calculations were performed using the GGA scheme. An 
expansion of Aa/a = 10% was used in the calculation (for Aa/a =0 see Sec lIII B"2l and also 
Ref.— ). This particular value was chosen as it is the case where the Fermi energy lies just 
in the middle of the gap of the minority DOS (see Fig0Ji,b) and thus the material would 
securely be in a HMF state. 

There are no minority states at ep confirming the compound to be in a HMF state. The 
high density below ep is dominated by d states being located at Co and Fe sites. The small 
density of states near €f in the majority DOS emerges from the strongly dispersing majority 
bands crossing the Fermi energy. From the spin resolved bands, it is seen that the majority 
bands cross or touch the Fermi energy (e^) in rather all directions of high symmetry. On 
the other hand, the minority bands exhibit a clear gap around ep. For Co2FeSi, the width of 
the gap is given by the energies of the highest occupied band at the T-point and the lowest 
unoccupied band at the T or X-point. The smaller value is found between T and X, thus it 
is an indirect gap. The conclusion drawn from the displayed electronic structure is that the 
states around the Fermi energy are strongly spin polarized and the system is according to the 
calculations indeed a half- metallic ferromagnet, at least for an enlarged lattice parameter. 
The gap in the minority DOS has the compelling result that the magnetic moment is integer 
with a value of QfiB, as expected from the experiment. 

The change of the electronic structure with the lattice parameter is clear from two basic 
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points of view. Firstly, from a free electron like band structure one expects a change of the 
band gradients and thus simply a shift of due to the accompanied changes in the DOS. 
Secondly, from a tight binding like approach one expects a difference in both, overlap and 
hopping integrals. In particular, they will be smaller for larger lattice parameter. With a 
decrease of a, the interaction between the atoms becomes stronger and the higher overlap 
results in a stronger de-localization of the electrons. Obviously the overlap is too large in the 
calculation using the experimental lattice parameter of Co2FeSi. On the other hand one may 
make use for nearly HMF compounds that small deviations of the magnetic moment from 
integer values can be compensated by small variations of the lattice parameter supposed 
the LSDA (GGA) calculations result in the correct electronic structure. Unfortunately, the 
opposite is also true. If the calculated Fermi energy is not just in the middle of the minority 
gap, then a small deviation of the lattice parameter may destroy the HMF character. 

The major conclusion that can be drawn from the needed increase of the lattice parameter 
of Co2FeSi is that bare LSDA (GGA) calculations are not sufficient to explain its electronic 
and geometric structure at the same time and thus may fail for similar compounds, too. 
Moreover, it is suggested that electron-electron correlation might play an important role for 
opening the gap in the minority states and to get the magnetic moment with an integer 
value. The verification of the experimental results has to be fulfilled independent of the 
material for both, lattice parameter and magnetic moment. In the following section it will 
be analyzed whether this goal can be reached possibly by inclusion of correlation. 

B. Electron correlation 

The relative importance of itinerant versus localized properties of d electrons in metal 
alloys was already discussed by van Vleck^ and Slate r 2 ^ . Heusler compounds are gener- 
ally thought as systems exhibiting localized magnetic moments as was first mentioned by 
Pauling for C^MnAl 26 . In particular, the magnetic moments of the C02YZ half-metallic 
ferromagnets follow strictly the localized part of the well known Slater-Pauling curved. It is 
clear that the d electrons are de-localized in metals. Thus, the question about correlation in 
transition metal based compounds is to what extent the on-site Coulomb interaction between 
d electrons is preserved such that important atomic properties like Hund's rule correlation 
are significant and determine the magnetic properties at least partially^. 
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The idea from the above discussion is, after the structural optimization failed to explain 
the magnetic moment of Co2FeSi, that inclusion of electron-electron correlation may be 
necessary in order to respect a partial localization of the d electrons in a better way. There- 
fore, the LDA+U method was used to recalculate the electronic structure of Co2MnSi and 
Co 2 FeSi. 

The LDA+?7 scheme is designed to model localized states when on-site Coulomb inter- 
action becomes important. It gives a self energy correction to localized states embedded in 
de-localized states. The energy U of the Coulomb interaction is rather large in free atoms 
(17eV to 27eV in 3d transition metals, see FigEJ), while screening results in solids in a much 
smaller values 4 ^ (for example 4.5eV in bcc Fe). In the present work, the value of the effective 
Coulomb-exchange interaction U e ff = U — J was varied in order to reproduce the measured 
magnetic moment. It turned out that U e fj from 2.5eV to 5.0eV for Co and simultaneously 
2.4eV to 4.8eV for Fe result in a magnetic moment of 6/i#. These values for U e ff corre- 
spond to about 7 . . . 20% of the free atom values as calculated from the corresponding Slater 
integrals by means of Cowan's program 4 ^ 

The average on-site Coulomb energy of pairs (ij) of equivalent electrons (i = j) is calcu- 
lated from the Slater integrals to be: 

En = F ou - -rj^-—J2 ] F *>,u (!) 

4l i + 1 ^ J 

The sum runs over all k > resulting in non-zero 3j-symbols. The 3j-symbols vanish for 
odd k and for d electrons one has k max — 2li — 4 such that EqJT] results in: 

2 2 

E d d = F - —F 2 - — F 4 (2) 

UO DO 

where the on-site Coulomb integral is U = F and all integrals are calculated for / = 2. The 
exchange integral is given by: 

J=^(F 2 + F 4 ) (3) 

Equation |H1 as calculated from Cowans's work 43 differs from the one reported by Anisimov 
et giving a factor 1/14 instead of 2/63, for unknown reasons*^. 

The atomic values of U = F , J, and U e ff = E d d as calculated for the 3d transition 
metals are displayed in FigJSJ In the following we will denote the U e ff values by U x where 
the subscript x stands for the values (in %) relative to the atomic ones (neutral atoms). 
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FIG. 5: Atomic Coulomb-exchange parameter. 

The values were calculated for the 3d transition metals using the neutral as well as the most 
common ionic configurations. Note that the values for J in (b) are shown on an enlarged scale 
below the break. 

It should be mentioned that the effect of screening in solids will mainly influence U rather 
than the exchange integrals J. On the other hand, the values of J depend also on the 
ionicity (see FigEb) that will be different in solids and dependent on the compound under 
investigation. From a practical point of view, the use of U x relative to the atomic values 
allows an easier comparison between different systems compared to the use of separate values 
for all quantities and atoms. 

First, a structural optimization was performed for C^FeSi using different values of U e ff 
in the calculations (see also Sec Jlll A|) . From Fig|T]it is seen that adding U e ff results in a less 
pronounced change of the £'(a)-dependence and shifts the energy minimum slightly to larger 
values of a. The minima are found at 5.72Aand 5.75Afor Uiq and U15, respectively. These 
values correspond to Aa/a exp of about +1.6% and thus are slightly higher compared to the 
experimental value. U e ff has the effect that the total energies at the experimental lattice 
parameter are about 3eV lower for {7 15 compared to Uio. At the same lattice parameter, the 
total energy of the LSDA calculation is about 6 eV higher compared to Uiq. It is worthwhile 
to note that the total energy taken from the LDA+/7 calculation does not correspond to 
the LDA ground state energy anymore and thus may lead to wrong conclusions about the 
structure. However, E(a, U e ff) will still give an idea about the optimal structure. 
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TABLE I: Magnetic moments of Co 2 FeSi. 

Given are the values calculated for a = 5.64Ausing different calculation schemes (see text). All 
values are given in \xb- Total moments {vo-tot) are given per unit cell and site resolved values (spin 
moment m s , orbital moment rm) are per atom. (U was set to 15% of the atomic values. +SO 
assigns calculations with spin-orbit interaction included. ) 





Co 




Fe 




Co 2 FeSi 




m s 


rm 


m s 


rm 


mot 


LSDA (vBH) 


1.31 




2.72 




5.29 


LSDA (VWN) 


1.40 




2.87 




5.59 


GGA (PBE) 


1.39 




2.85 




5.56 


GGA + SO 


1.38 


0.04 


2.83 


0.06 


5.56 


LDA+J7 


1.53 




3.25 




6.0 


LBA+U + SO 


1.56 


0.08 


3.24 


0.07 


6.0 



1. Magnetic moment and minority gap in LDA+U 

In the following the dependence of the magnetic moments on the type of the exchange- 
correlation functional is analyzed and compared to the experimental values. The results of 
the calculations for Co2FeSi using different approximations for the potential as well as the 
parametrization of the exchange-correlation part are summarized in Tab|lJ 

Tabel|l]gives the site resolved moments at Co and Fe sites and the total magnetic moment. 
The induced moment at the Si sites (not given in TabUwas aligned in all cases anti-parallel to 
that at the transition metal sites. Likewise, the magnetic moment located in the interstitial 
is omitted in Tab|lJ however, both quantities have to be respected in order to find the correct 
total moment. 

The lowest value for the total magnetic moment was found if using the v. Barth-Hedin 
(vBH)22 parametrization with a discrepancy of -0.7/iB compared to the experiment. The 
results from Vosko-Wilk-Nussaii 1 ^ parametrization and GGA— are very similar but still too 
low compared to the experiment. It is seen that both site specific moments are too low to 
reach the experimental value that may need values of about 1.5/x.b and 3/i# per atom at 
the Co and Fe sites, respectively. From this point of view, the LSDA or GGA approaches 
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are not sufficient to explain the magnetic structure of Co2FeSi, independent of the type of 
parametrization. 

The correct magnetic moment at the experimental lattice parameter was only found if 
using the +U functional. The use of the LDA+[7 scheme improves the total magnetic 
moment considerably (C/15 was used for the calculation in Tab|lJ). The ratio of the magnetic 
moments of Co and Fe was measured to be mpe/^Co = 2.2 at 300K in an induction field 
of 0.4T— . The ratio of 2.1 found in the LDA+U calculations agrees very well with this 
value. Spin-orbit interaction (SO) was additionally included in some of the calculations to 
analyze its influence on the total and partial magnetic moments. In the GGA calculations, 
it did not improve the total moment. The experiments delivered orbital to spin magnetic 
moment ratios (mi/m s ) of about 0.05 for Fe and 0.1 for Co. The LDA+[/+SO calculations 
revealed 0.02 for Fe and 0.05 for Co, being a factor of 2 smaller compared to the experiment. 
Overall, the agreement of both ratios imp e jmc and m//m s ) with the experiment is better in 
LDA+[7 compared to the pure LSDA or GGA parametrization of the exchange-correlation 
functional. This may give advice that U corrects, at least partially, the missing orbital 
dependence of the potential in LSDA. 

In the following the effect of the magnitude of the Coulomb-exchange interaction on the 
magnetic moments will be discussed. Figure El compares the dependence of the magnetic 
moment on U e ff for Q^FeSi and Q^MnSi. 

It is found for Co2FeSi that the total magnetic moment increases from 5.6/is to 6/i# while 
increasing U e ff from to U7.5 and stays at the integer value up to U20. Above that value 
the moment exhibits a further increase. The situation is different for Co2MnSi with 5/xb at 
U e ff = 0. The moment stays integer up to about t/7.5 and then increases with increasing 

Ueff- 

In Co2MnSi, the moment at the Mn sites increases nearly linearly with U e ff. At the same 
time, the moment decreases at the Co sites up to U±q counter-balancing the Mn moments. 
At higher U e ff, the Co moment increases and the HMF character is lost. The situation 
is different in Co2FeSi, both site specific moments show an overall slight increase with 
increasing U e ff. Moreover, they stay rather constant in the range of U7.5 to U20 such that 
the total moment stays integer in this region. In both materials, the moment at the Co 
sites reacts stronger on U e ff than the ones at Y sites for U e ff being above the region of 
integer total moments. The increase of the site specific moments is much less pronounced 
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a) Co 2 MnSi b) Co 2 FeSi 




effective Coulomb parameter £/ ff / U [%] 

FIG. 6: Dependence of the magnetic moments on U e ff. 

The dashed horizontal lines assign the experimental values of the magnetic moments for Co2MnSi 
(a) and Co2FeSi (b). The calculated values are connected by solid lines to guide the eye. 

in Cc^FeSi compared to Co2MnSi. This gives advice that the Fe compound is much more 
stable against correlation effects in comparison to the Mn compound, at least concerning 
the magnetic moments. 

Again, it is interesting to find out whether the integer magnetic moment of Co2FeSi is 
related to a real HMF gap in the minority band structure and to compare the behavior of 
the gap to the one of Co2MnSi. The extremal energies of the gap enveloping bands are 
shown for both compounds in Fig|7| 

In both cases the size of the gap increases with increasing U e ff. From FigEK it is found 
that Co2MnSi stays in an HMF state in the range up to about Us, that is the range of the 
integer magnetic moment. Taking larger values shifts ep to be outside the gap. In Co2FeSi 
(Figdo) , the minority gap involves ep from about U% up to L^o- That means, the integer 
value of the magnetic moment is related to the minority gap in both compounds and thus 
the direct consequence of the HMF state. 

In both compounds, the gap is completely destroyed at very large values of U e ff (> 8eV). 
The effect of the Coulomb-exchange interaction on single atoms (Co, Mn, or Fe) on the 
minority band gap was also investigated. It was found that the minority gap is destroyed in 
both materials if U e ff is added only to one of the 3d elements, as expected. However, small 
changes in only one of the U e ff values about the balanced value did not change the general 
behavior. 
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FIG. 7: Dependence of the minority gap on the effective Coulomb-exchange parameter. 

Shown are the extremal energies of the gap involving states for Co2MnSi (a) and Co2FeSi (b). The 

shaded areas assign the region of half-metallic ferromagnetism. Lines are drawn to guide the eye. 

2. Electronic structure in LDA+U 

In the following the influence of the correlation on the electronic structure will be discussed 
in more detail. Figure |H1 compares the spin resolved density of states for C^MnSi and 
Co2FeSi calculated in LSDA and LDA+£7 approximation. The low lying sp-bands emerging 
from T^g (majority) and Y§ g (minority) s states are located at (—8. . . — 10)eV. They are 
not shown here (compare FigEJ). These states are coupled to p states (r 5u 6u ) and mainly 
located at the Si atoms. 

In both materials, the majority d states spread from 8eV below to about 4eV above ep 
where the onset of the high lying, unoccupied s states becomes visible, s states are also 
present at below -4eV due to the s — d coupling of the T 5g states, p states are contributing 
to the density over the hole range of the d bands. They are mixed to s and d states due to 
the coupling of states with even and odd parity if going away from the T-point. The high 
density of states at -leV in Co2MnSi or -2eV in Co2FeSi emerges mainly from r 89 like states. 

The situation is similar in the minority density of states where additionally a gap exists. 
This gap splits the minority rf-states and the unoccupied part above ep is mainly of Ti g 



16 



a) Co,MnSi b) Co 2 FeSi 




energy E-£ f [eV] energy E-e p [eV] 

FIG. 8: Spin resolved densities of states for Co 2 MnSi and Co2FeSi. 

Gray lines indicate the DOS within LDA+C using f7i for Co2MnSi (a) and for Co2FeSi (b). 
Black lines indicate the results from LSDA. The upper and lower parts of the plots display the 
majority and minority states, respectively. (See text for the particular values of U at the different 
sites.) 

character. The minority gap is clearly different in both compounds and depends strongly 
on the applied effective Coulomb-exchange interaction as already discussed above. 

In the case of Co 2 FeSi, where the experiment clearly indicates the total magnetic moment 
to have an integer value, the LSDA results seems to be incorrect. On the other hand, it 
should be noted that the unoccupied d states are sharply peaked at €f, resulting in a DOS 
that would imply an inherently unstable system in usual, simple ferromagnetic materials. 
This is despite the spin polarization. The system is, however, stabilized by the gap in the 
minority states that fixes the number of occupied states not only in the minority but also - 
and more important - in the majority channel. In LDA, the paramagnetic e g (located at the 
Co sites) and t 2g (at Fe) states are rather narrow at ep suggesting that Coulomb correlation 
of the Hubbard U kind, being ignored within LSDA, will make Co 2 FeSi insulating in the 
minority states. Attempting to treat electron correlation, the electrons become stronger 
localized by applying the LDA+U method. One of the major issues here is that the LDA+£7 
type correlation is important for the electronic and magnetic structure of Co2FeSi. In more 
detail, the majority states at around -4 eV (arising from t 2g states in the paramagnetic case) 
are shifted by C/15 to -6eV below ejr. At the same time the unoccupied minority states seen 
in the LSDA calculation just above exhibit a larger shift of the paramagnetic e g states 
(to 2 eV above €f) compared to the paramagnetic t 2g states (to 0.5 eV above 6f)- The 
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occupied states of the minority DOS exhibit a less complicated behavior, they are mainly 
shifted towards the Fermi energy. 

For Co2FeSi as shown in FiglHt, U15 was chosen because for this value the Fermi energy 
lies just in the middle of the gap. It is clearly seen that the splitting between the occupied T 8g 
and the unoccupied T 7g states becomes larger. The minority T^g states exhibit an additional 
splitting after applying +U, whereas the unoccupied majority states stay rather unaffected. 
At the same time the gap in the minority states becomes considerably larger and the Fermi 
energy is clearly within this gap confirming that the compound is in a HMF state. 

In FigJEk the DOS of Co2MnSi is compared in LSDA and LDA+U approximation using 
Uiq being above the limit of the HMF state. The majority d states being mainly located 
at the Co and Mn sites are clearly seen below ep and possess a width of about 7eV. Near 
€f one finds mainly a low density of majority states emerging from strongly dispersing d 
bands. In LSDA, there are hardly any minority states at ep confirming the compound to be 
a half-metallic ferromagnet. Using LDA+U, the d densities are shifting away from ep similar 
to Co2FeSi. However, the occupied part of the minority DOS also changes if applying U e ff\ 
it shifts closer to the Fermi energy and crosses it for high values of U e ff. Once it touches ep 
the maximum value of is reached that still results in the HMF state what is here the case 
for about t/7.5. 

Figure M shows the changes of the band structure of Co 2 FeSi if applying the LDA+U 
method. For easier comparison, details of the band structures are shown along T — X that 
is the A-direction of the paramagnetic state. The A-direction is perpendicular to the C02 
(lOO)-planes. As was shown earlie^, just the A-direction plays the important role for the 
understanding of the HMF character and magnetic properties of Heusler compounds as was 
also pointed out by Ogiit and Rabe^. 

The developing of the band structure with increasing U e ff is shown in FigEl using values 
of 1.26eV, 2.52eV, and 3.78eV at the Co sites and simultaneously 1.2eV, 2.38eV, and 3.56eV 
at the Fe sites (U5, £Ao, and U15). The upper panels (a-d) display majority and the lower 
panels (e-h) the minority band structures. It is seen that only slight changes appear in the 
majority band structure, the width of the visible bands at T stays rather unaffected and the 
main changes appear close to X. 

More interesting is the behavior of the minority bands as those determine the HMF 
character of the compound. One finds that the energies of the occupied states at V stay 
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FIG. 9: Band structure of Cc^FeSi with variation of the Coulomb-exchange interaction. The gray 
lines assign the upper and lower bands defining the gap at the X-point in the minority states. (See 
text for the particular values of U% at the different sites.) 

nearly the same. The shapes of the bands close to T are similar, too. The situation is 
different at X where the unoccupied states are shifted away from €f resulting in a gap that 
increases with U e ff. It is clearly seen that it is an indirect gap. At C/ 10 the Fermi energy 
falls within the gap and the HMF state is reached. The minority gap is larger at T or at 
other points of high symmetry (not shown here) compared to X. Therefore, it is obvious 
that the states at X determine the HMF character of the compound. 

The situation is clear at this stage: the gap in the minority states stays only within a 
certain range of the effective Coulomb-exchange interaction. If one goes away from that 
limit the Fermi energy falls not longer inside of the gap and the material loses its HMF 
character. It is also obvious that the mechanism leading to the half-metallic ferromagnetism 
in Co 2 FeSi serves to destroy it in Co 2 MnSi. The worst case would appear for both materials 
at about 7... 8% of the atomic values of C/ e // where LDA+£7 predicts still values very 
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close to the measured magnetic moments but also is the borderline for the loss of the HMF 
character. From this point of view, nearly integer magnetic moments alone do not verify the 
half-metallic ferromagnetism and one may have to search for alternative materials. 

Previous band structure calculations using LSDA predicted Q^MnSi to be a HMF at the 
experimental lattice paramete r 24 ! 47 . Here, LDA+U calculations were performed in order to 
consider the influence of on-site electron correlation in this compound. It was found that 
the gap in the minority DOS stays up to U e ff = 2.3eV for Mn and 2.5eV for Co. For larger 
values of U e ff, that means for stronger correlation, the system loses its HMF character, by 
shifting €f outside the minority gap. C^MnSi is expected to exhibit 100% spin polarization 
at the Fermi energy. In reality this was not verified up to now. Wang et a& found only about 
10% spin polarization at room temperature from a spin-resolved photo emission experiment. 
The spectra shown in Ref.— give no advice on the existence of a surface state that may serve 
to lower the spin polarization at ep. Here, it was found that an already moderate correlation 
will destroy the gap. The spin polarization at ep amounts to 75% if applying U\o, a value 
that is compatible to about 55% found in Ref.!^ 7 -. From that point of view, on-site correlation 
might be one more reason why no complete spin polarization was found in this compound 
(for others see^). 

Many effects have been considered to explain the missing experimental proof of the com- 
plete spin polarization in predicted half-metallic ferromagnets^. For sure, each of these 
effects may serve by itself to reduce the spin polarization at the Fermi energy. On the 
other hand, none of these effects is able to explain a difference of 15% between the ob- 
served and calculated magnetic moments like for the case of Co2FeSi. The appearance of 
surface or interface state a l5 ° may decrease the spin-polarization in certain layers but will 
not influence drastically the measurements of the magnetization of bulk samples. Magnon 
excitation a 51 i 52 i 53 i 54 i 55 will lead to a decrease of the magnetic moment with temperature 
rather than in an increase and thus can easily been disregarded here. Spin-orbit interaction 
couples, indeed, spin- up and spin-down states^ but does not enhance the magnetic moment, 
as was shown above. From this point of view, correlation might be one of the main causes of 
the too low spin polarization being observed in experiments. In the presented work, a static 
correction was used in the LDA+£7 approximation. Dynamic correlation effects may be 
respected by means of the LDA+DMFT approximation using dynamical mean field theory 
(DMFT). Using this method for NiMnSb, it was shown by Chioncel et af 1 that non quasi- 
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particle effects may also serve to change the properties of the gap in the minority states. 
The use of bulk sensitive angle-, energy-, and spin- resolved photo emission becomes highly 
desirable in order to explain the details of the band structure in the most complete experi- 
mental way and to proof the existence of correlation in predicted half-metallic ferromagnets 
on the one hand and the true existence of half-metallic ferromagnetism on the other. 

IV. SUMMARY AND CONCLUSIONS 

Structural parameters, magnetic moments and electronic properties of the Heusler com- 
pounds Co 2 MnSi and Co 2 FeSi were presented. The obtained, optimized lattice parameter 
leads for Co 2 FeSi to a too small magnetic moment, using LSDA or GGA calculations. The 
measured value of the magnetic moment was only found for a lattice parameter enlarged 
by about 10% compared to the experiment. It was also shown that the magnetic moment 
of Co 2 MnSi is more stable against structural changes in the LSDA or GGA approaches. 
At present, it seems that simple LSDA and GGA methods are not sufficient to explain the 
electronic structure of Heusler compounds completely, at least they fail for Co 2 FeSi. 

In the step beyond LSDA, it was shown for the case of Co 2 FeSi that electron correlation 
is able to explain the experimental magnetic moment at the experimentally observed lattice 
parameter. It was found that the LDA+£7 scheme reproduces satisfactorily the experimental 
observations. At moderate effective Coulomb-exchange correlation energies of about 2.5eV to 
4.5eV, the LDA+[7 calculations agreed very well with the measured total and site specific 
magnetic moments. At the same time the compound was predicted to be clearly a half- 
metallic ferromagnet. 

Finally, it was checked how the inclusion of correlation will change the calculated prop- 
erties of similar Heusler compounds. Co 2 MnSi is a compound being close to Co 2 FeSi in 
the series of known Heusler compounds based on Co 2 . It is obvious that on-site correlation 
will also play an important role in the Mn compound, if it does in the Fe compound. It 
was found that a small correlation energy of only about 2.5eV destroys the half-metallic 
ferromagnetism what may explain why a complete spin polarization was not observed in 
this compound up to now. 

In conclusion, it is suggested that it will be necessary to include electron-electron corre- 
lation beyond LSDA and GGA in the theoretical description of all potential half-metallic 



21 



ferromagnets and in particular for Heusler compounds. 
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